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Abstract 

Collocation implementation of the Kleiser-Schumann’s method in geometries with two non- 
periodic directions is shown to suffer from three spurious modes - line, column and checkerboard 
- contaminating the computed pressure field. The comer spurious modes are also present but they 
do not affect evaluation of pressure related quantities. A simple methodology in the inversion of the 
influence matrix will efficiently filter out these spurious modes. 
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1. Introduction 


Most existing numerical solutions of incompressible Navier-Stokes equations in three di- 
mensions employ the primitive-variable formulation. Here, the velocity field and pressure can not 
be approximated independently and must satisfy a compatibility condition. Approximating both ve- 
locity and pressure by polynomials of the same degree will result in some spurious modes for the 
pressure [1]. In a spectral collocation implementation, these spurious modes can be characterized 
as the spurious components of the pressure field whose discrete gradient at the interior collocation 
points, where discretized momentum equations are satisfied, is zero. Such pressure components 
have no effect on velocity and are therefore left uncontrolled by the discretized governing equations. 

Existence of spurious pressure modes in a spectral simulation was first pointed out by Mor- 
choisne [2]. In coupled spectral implementations where the continuity equation is discretized direct- 
ly, theoretical analysis of the spurious modes is on a firm footing [1]. When a non-staggered grid 
is employed it has been shown that a fully periodic problem has no spurious modes, discounting the 
arbitrariness of the mean value of pressure in an incompressible flow. In a flow with one non-period- 
ic direction pressure has one spurious mode, and problems with two non-periodic directions have 
seven ( l -line , \-column, 1 - checkerboard and 4 -corner) spurious modes [1,3]. These results have 
been shown to apply to the following spectral implementations: Galerkin Legendre, Legendre tau, 
Galerkin Chebyshev, Chebyshev tau and Chebyshev collocation with Gauss-Lobatto points. In a 
collocation implementation, the spurious modes can be avoided with an appropriately staggered 
mesh for velocity and pressure. Montigny-Rannou and Morchoisne [4] have recently described an 
algorithm for the two non-periodic problem on a half-staggered grid which involves only one 
0 checkerboard ) spurious mode. The spurious modes can be completely avoided in a two non-period- 
ic problem with a fully-staggered mesh [5, 6]. 

A second class of numerical algorithms avoids direct solution of the continuity equation by 
solving a Poisson equation for pressure [ 1 ] . The time-split implementation [7] is the simplest of them 
all but it suffers from non-zero boundary divergence, non-zero slip velocity and time-splitting errors. 
Kleiser- Schumann’s influence matrix method can be used to enforce zero boundary divergence and 
no-slip condition [8,9]. A partial implementation of the influence matrix is simpler, but will result 
in small but non-zero interior divergence [10]. A full implementation of the Kleiser-Schumann’s 
method with tau (or collocation ) correction enforces interior zero divergence as well. Here we iden- 
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tify the spurious modes in the collocation implementation of the Kleiser-Schumann method [8, 9, 
1 1] with tau correction for a problem with two non-periodic directions to be the line, column, check- 
erboard and comer modes. A simple correction procedure which will automatically filter these spu- 
rious modes is obtained. This correction procedure is applied in the simulation of a turbulent square 
duct flow [11] and is found to be very effective in eliminating spurious modes. 

2. Spurious Modes for Kleiser-Schumann’s Method 

The full implementation of the Kleiser-Schumann’s method with collocation correction can 
be given by the following equations and boundary conditions: 


V 2 p = 

- V ■ (NL) + V • B 

in 
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CL»|3f 
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- -Vp + ^V J V + B 

in 

Q,dQ 


II 
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B = 0 
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where NL is the nonlinear term in the Navier- Stokes equation and Vf, is the velocity boundary condi- 
tion. The above equations and boundary conditions are in their discretized form, therefore the sym- 
bols V,V • and V 2 represent discrete gradient, divergence and Laplacian operators, and Q and dQ 
represent interior and boundary collocation points. Since the momentum equation is satisfied only 
in the interior, B represents the boundary momentum residual and is nonzero only on the boundary. 
Combining equation (la) and the discrete divergence of equation (lb) one obtains 

' Y. = V 2 (V • V)\nQ, provided the discrete divergence of the discrete gradient operator is 
ot Re 

identically equal to the discrete Laplacian operator. The above equation for the velocity divergence 
along with the boundary condition given by equation Id, results in a divergence free velocity field 
both in the interior and on the boundary. 

By definition, each spurious mode is a valid solution to the discretized governing equations 
and appropriate boundary conditions. The spurious modes have a non-zero contribution to pressure 
but have no effect on velocity, therefore V sp = 0, where subscript “sp” stands for the spurious 
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mode. The spurious pressure components therefore satisfy 

V 2 Psp = VB sp in Q 

0 = - V Psp + B sp in Q,dQ 


B sp = 0 m Q 

where B sp is the corresponding spurious boundary momentum residual. Any non-trivial solution to 
the above linear equations represents a spurious mode, which when added to the true solution will 
still satisfy the discretized Navier-Stokes equations and boundary conditions. 

Eight solutions to the above equations can be identified. First of which is the non-spurious 
solution, pip = constant and B sp = 0, which indicates that pressure is evaluated only up to an arbitrary 
additive constant in incompressible flows. The first two spurious modes are the line and colutnn 
modes 


P sp 


T Nx (x) i 


Psp = T Ny (y) , 



1 ) Nx (Nx) 2 


lA/Vy ) 2 


at x = ± 1 
otherwise 

at y = ± 1 
otherwise 


B 2sp sp 0 

B w = B }sp = 0 


( 3 ) 


Here (Nx+ 1) and (Ny+\) are the number of points along the non-periodic Chebyshev directions and 


the third direction is at most periodic. B } , B 2 and B 3 are the three components of the boundary mo- 
mentum residual and only the normal component is nonzero. T^ x Tjyy are Chebyshev polyno- 
mials of the highest degree along x and y. The third spurious mode is the checkerboard mode, with 
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These three spurious modes have no variation in the periodic z direction and therefore contaminate 
only the zero* mode along the z direction. The other four spurious modes are the comer modes and 
each of them can have arbitrary variation along the z direction. For example, let be the arbi- 
trary variation along the x=l, y= 1 corner. Then the comer mode corresponding to this comer can 


= o 


now be written as 
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where D/^ x (x) is the discrete derivative of the polynomial which collocates to zero at all points except 
at je=l. The comer modes for the other three comers can be written similarly. The comer spurious 
modes simply reflect the fact that in a collocation implementation the pressure along the four corner- 
lines never enter into the computation and therefore their values remain unspecified. This arbitrari- 
ness associated with the corner pressure values is relatively innocuous, since they are not required 
in computing pressure related quantities of interest. On the other hand, the other three spurious 
modes are buried in the numerically computed pressure and need to be filtered out. 


3. Filtering Procedure 


Implementation of the Kleiser-Schumann’s method involves the construction of an influence 
matrix. The solution of the pressure Poisson equation (equation la) requires the knowledge of pres- 
sure boundary conditions ( pb ) and boundary momentum residuals (Bb) at the (2N x +2N y -4) points, 
excluding the corner points. Note that in the three, line, column and checkerboard modes only the 
normal component of the boundary momentum residual is nonzero. In essence, a total of 
(4N x +4N y -%) unknown boundary pressure and normal momentum residuals are required in evaluat- 
ing the pressure field. These unknowns are evaluated by requiring that continuity and normal mo- 
mentum (with the residual) are satisfied at the boundary points excluding the corner points. This pro- 
vides the necessary (4A^+4A^-8) linear equations for the unknown quantities. These equations can 
be cast into the following matrix form, A x = R, where A is the influence matrix, i is the unknown 

vector of boundary pressure and normal momentum residual and R is the right hand side. In a three- 
dimensional problem, a Fourier transform along the periodic z direction will result in one influence 
matrix for each Fourier wavenumber k z . Invertibility of the influence matrix is closely related to the 
presence of spurious modes. In particular the influence matrix corresponding to k z =0 suffers from 
the constant, line, column and checkerboard modes. This influence matrix therefore has four zero 
eigenvalues. In the numerical computation of a turbulent square duct flow it was observed that only 
for 0 were four eigenvalues of the influence matrix of the order 10 -12 , while the lowest of all other 
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eigenvalues were around 10 -6 . The eigenvectors corresponding to these eigenvalues are the corre- 
sponding spurious boundary pressure and normal component of the momentum residual. 

As suggested by Tuckerman [9], the non-invertibility of the influence matrix can be easily 
overcome by constructing a related matrix A' = M A' M *, where M and M 1 are the eigenvec- 
tor matrix of the original influence matrix and its inverse, respectively, and A’ is a diagonal matrix 
with the eigenvalues along the diagonal and with the zero eigenvalues replaced by some non-zero 

constant. The influence matrix now becomes invertible, i.e., 1 = M 1/A' M \ and the result- 

ing pressure and boundary momentum residuals yield a divergence-free flow field independent of 
the constant that replaces the zero eigenvalue. Let the p lh eigenvalue of the original influence matrix 
be zero and be replaced by a constant c p . Let the corresponding p lh eigenvector be M ip , which is a 
vector of boundary pressures and momentum residuals corresponding to a linear combination of the 

Ad ipbp 

spurious modes. Contribution of this p lh mode to the unknown vector*; is then — - — , where b p is 

c p 

the projection of the right hand side along the eigenvectors, given by M ~ 1 Rj. One simple way to 
filter the four spurious modes will then be to set the arbitrary constant c p to be infinity. In other words, 
in the evaluation of A^\ one over the zero eigenvalue is replaced by zero. This filtering procedure 

was implemented in the computation of turbulent flow in a square duct [11] and it essentially re- 
moved all high frequency checkerboard type oscillations. It should be pointed out that such oscilla- 
tions can be observed in the k z ^() modes as well, especially at high Reynolds numbers when resolu- 
tion is only marginal. These oscillations are not due to the spurious modes, but are controlled energy 
in the high wavenumber modes. 

4. Spurious Modes in Partial and Time-Split Methods 

Spurious modes for the partial implementation of the Kleiser-Schumann’s method without 
the collocation correction can be analyzed in similar fashion. The spurious modes satisfy the follow- 
ing equations. 
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V 2 Psp — o 
0 = — Vpsp + B sp 


in Q 


(6a-6c) 


in Q,dQ 


B S p - 0 in Q 

It can be easily seen that the only admissible solution for the above set of equations is the non-spuri- 
ous constant mode and there are no spurious modes present. This is confirmed in the numerical simu- 
lation by observing that the influence matrix for the k z =0 mode has only one zero eigenvalue and 
the corresponding constant mean pressure can be set to zero by replacing the zero eigenvalue by in- 
finity. 

Spurious modes in the time split implementation will depend on the exact boundary condi- 
tions employed for the intermediate star-level velocities (V*) and the pressure Poisson equation. Fol- 
lowing Streett and Hussaini [7], if we employ [(2 Vp(r) — Vp(t — 6tj) • r] as the boundary condition 

for the tangential components of the intermediate star-level velocity, zero penetration for the normal 
velocity component and zero Neumann boundary condition for the pressure, then we have the fol- 
lowing equations satisfied by the spurious modes: 


V* p = At Vp sp 



y%'v = Vpsp • v = o 


in Q 


in Q 
in bQ 


(7a-7d) 


V* S p • r = ( 2Vp sp (t ) - Vpsp{t - <5r)) • r in dQ 


where rj and r are direction normal and tangential to the boundary. The analysis of the spurious 
modes is more complicated and also depends on the initial tangential pressure gradients on the 
boundary. With careful choice of initial conditions, the no penetration and pure Neumann boundary 
conditions will guarantee no spurious components. 


5. Conclusion 


Collocation implementation of the Kleiser-Schumann’s method in geometries with two non- 
periodic directions have three spurious modes - line, column and checkerboard- contaminating the 
computed pressure field. The corner spurious modes are also present but they do not affect evalua- 
tion of pressure related quantities. The three spurious modes can be easily filtered out by replacing 
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the zero eigenvalues of the influence matrix with infinity before solving for the unknown boundary 
pressure and momentum residuals. Partial implementation of the Kleiser-schumann’s method with- 
out collocation correction admits no spurious modes. Spurious modes can also be avoided in time- 
split implementations. 
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